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Predicted masses of J p — + and 1 + states of D s j and J) by a potential model proposed some 
•"T' ■ time ago by two of us (T.M. and T.M.), in which the Hamiltonian and wave functions are expanded 

in 1/mQ with tuq heavy quark mass respecting heavy quark symmetry, have recently been confirmed 
f*"^ ' by BaBar and Belle experiments within one percent accuracy. 

In this Letter, decay modes of + and 1 + states of B and B s mesons are discussed using the 
predicted masses of this model. 
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"sj- ■ I. INTRODUCTION 

o 

The BaBar's discovery [1] of D s j(2317) and D s j(2460) of the system c and s quarks has inspired theorists [2, 3] 
to explain these states in some way because a well-known potential model [4, 5] failed to expound masses of these 
states. This discovery has been soon confirmed by CLEO and Belle [6]. These D s j states are narrow because isospin 
violation and DK/D*K threshold prohibit them to decay into lower states D s (1968) and D*(2112). This is another 
• reason why theorists are interested in these states and are trying to explain their branching ratios, too. 

Soon after the discovery, another set of heavy mesons, £>o°(2308) and D 1 °(2427) of c and u/d quarks which have 
the same quantum numbers J p = + and 1+ as D s j, has been discovered by Belle [7]. Again masses of these states 
cannot be fitted with those of a potential model [4, 5]. These states Dq° and have no restriction like D s j and 
■ hence their decay (like Dn and D*ir) width becomes broad. 

To explain masses of D s j, Bardeen, Eichten, Hill, and others [2, 3] proposed an interesting idea of an effective 
Lagrangian with chiral symmetries of light quarks and heavy quark symmetry. The heavy meson states with the 
total angular momentum J = and J = 1 with j q = 1/2, which is the total angular momentum of the light quark 
degrees of freedom, make the parity doublets (0~,0 + ) and (l - , 1 + ), respectively, and the members in each doublet 
degenerate in a limit of chiral symmetry. Furthermore, the two states (CP, 1 _ ) degenerate in a limit of heavy quark 
symmetry, as well as (0 + , 1 + ) do. These doublets, (0~, 1 _ ) and (0 + , 1 + ), are called the heavy spin multiplets. Though 
their model can reproduce well the mass difference between members in each parity doublet by using the modified 
Goldberger-Treiman relation, unfortunately it cannot calculate masses themselves. This is why some people have 
proposed a modified potential model by using a DK bound state [8]. 

Several years ago, Godfrey, Isgur, and Kokoski proposed a relativised potential model [4, 5], in which they included 
\/ p 2 + m 2 as a kinetic term and spin-dependent interaction terms where linear-rising confining as well as short-range 
Coulomb potentials are taken into account. This model, however, could not reproduce the masses of D s j(2317) and 
D s j(2460). Even though they have taken an infinite heavy quark mass limit in the second paper [5], they could 
not reproduce D s j(l + ) mass and their calculated masses exceeded DK/D*K threshold. Thus, there appear some 
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TABLE I: Comparison of higher D and D s meson masses (units in MeV) 



J p 


D(0+) 


D(l+) 


D s (0+) 


D.(l+) 


observed 


2308 


2427 


2317 


2457 


predicted 


2304 


2449 


2339 


2487 



TABLE II: Higher B and B a meson masses taken from [9] (units in MeV) 



J p 


B(0+) 


B(l+) 


Bs(l-) 


B 3 (0+) 


Bs(l + ) 


predicted 


5697 


5740 


5440 


5716 


5760 



discussions [8] that a potential model is not appropriate to describe these states. 

The above potential model does not completely and consistently respect heavy quark symmetry and does not treat 
quarks as Dirac particles. Some time ago, two of the authors (T.M. and T.M.) [9] proposed a new bound state 
equation for atomlike mesons, i.e., heavy mesons composed of a heavy quark and a light antiquark. In this model 
compared with [4, 5], quarks are treated as four-spinor particles from the beginning and both the Hamiltonian and 
wave functions are expanded in 1/ itiq so that our model treats quarks as relativistic as possible and consistently takes 
into account heavy quark symmetry within a potential model. Our predicted masses for these states, D s j and Dq° 
and D] , are in good agreement within one percent accuracy when calculated at the first order in 1/toq. See Table I 
in this letter and Tables III and IV in [9]. 

In this Letter, having confirmed that our model has well succeeded in predicting masses of recently discovered heavy 
mesons, we predict masses of + and 1+ of B and + , 1~ and 1+ of B s mesons by citing the evaluated values in [9]. 
Here we also discuss their decay modes; whether these mesons violate isospin symmetry or not and whether mass 
difference between + (1 + ) and 0~ (1 _ ) is less than the BK / B* K-m&ss threshold or not. We expect these higher 
states of B and B s mesons can be detected in CDF/LHC experiments. 

II. HIGHER B AND B a MESON MASSES AND THEIR DECAY MODES 

Our prediction in [9] of B and B s meson masses in the first order of 1/ toq corrections is given in Table II. Considering 
the fact that D/D s masses have been well predicted by our model within one percent accuracy as seen in Table I, we 
expect masses of these B and B s mesons are similarly within one percent accuracy to be observed. Furthermore when 
looking at Table I carefully, we expect that our predicted masses of B/B s mesons agree with experiments better than 
those of D s since m u ^/mb < m u ^/m c < m s /mb < m s /m c <C 1 (m Ui <j = 10 MeV, m s — 94.72 MeV, m c = 1457 MeV, 
and mb = 4870 MeV are adopted in our paper [9]). 

Let us discuss decay modes of + and 1 + states of B and B s mesons taking into account these meson masses given 
by Table II. 



(1) 


Bs(l-) 


^B s ((T)+7 




(2) 


B(0+) - 


-> B(O-) +7T 


with broad decay width 


(3) 


5(1+)- 


-> B(l-)+7T 


with broad decay width 


(4) 


B s (0+) 


-^B s (0-)+n 


with narrow decay width 


(5) 


B s (l + ) 


-^B s (l-)+n 


with narrow decay width 



Comments are given as follows. 



(1) Decay B s (l ) — > B s (0 ) + 7 is similar to B(l ) — > B(0 ) + 7, i.e. B* — > B + 7, and is dominant decay mode. 

(2) Decay width of B(0+) -» B(0~) + tt is as broad as a few hundred MeV like D{0+) -> D(0~) + it [7] because this 
is strong decay and is not prohibited by isospin invariance since I(B(0~)) = I(B(0 + )) = 1/2 while I(tt) = 1, 
where I(X) is isospin of a particle X. 

(3) Decay width of B(l+) -> B{1~) + n is also as broad as a few hundred MeV like D(l+) -> D(l~) +ir [7] because 
this is strong decay and is not prohibited by isospin invariance since I(B(1~)) = I(B(1 + )) = 1/2 while I(ir) = 1. 
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(4) Decay width of B s (0 + ) — > _B s (0 _ )+7r is expected very narrow, a few MeV like the decay of D s (0 + ) — ► D s (0 _ )+7r 
since this decay mode is prohibited by isospin invariance due to I(B s (0 + )) = I(B S (Q~)) = while I(ir) = 1 and 
the predicted mass of B s (0 + ) is below BK threshold. 

(5) Decay width of B s (\ + ) — ► B s (l~) + n is also expected very narrow, a few MeV like the decay of D s (l + ) — > 
-D s (l~) + 7T since this decay mode is prohibited by isospin invariance due to I(B S (1 + )) = I(B S (1~)) = while 
I(tt) = 1 and the predicted mass of B s (l + ) is below B*K threshold. 

Here B s (0~) = 5,(5370), D a (0~) = D±(1968), D s (l~) = D s *±(2112), D s (0+) = £ sJ (2317), and D s (l+) = D sJ (2460) 
[10]. We expect these higher states of B and B s mesons can be detected in CDF/LHC experiments by looking at 
their decay modes. 



III. CHIRAL LIMIT OF H 

As mentioned above, the model by Bardeen and Hill [3] has chiral symmetry in the chiral limit of light quark mass, 
m q — > 0. Our model has also chiral symmetry. Here let us see how it is realized. 

In our model the lowest order mass of the Qq bound state is given by tuq + Eq after solving the following Schrodingcr 
equation [9, 11], 

H ®Vo=£oV>o. H Q = a q -p + p q (m q + S(r)) + V(r), (1) 

where H is the lowest order Hamiltonian, I expresses all the quantum numbers, tp^ is the lowest solution to Eq. (1), 
and with a notation (g> one should understand that gamma matrices for a light antiquark be multiplied from left with 
the wave function while those for a heavy quark from right. Here S(r) is a confining scalar potential and V(r) is a 
Coulombic vector potential at short distances. Both potentials have dependency only on r, relative distance between 
Q and q. Quantities with a subscript q mean those for a light antiquark. Using the spherical polynomials Y™ and 
the vector spherical harmonics defined by 

we can decompose the wave function ipQ into radial and angular parts as follows [9, 11]; 

k 



where the angular part yf^ is given by the linear combination of Y™ and a ■ Yj m (A=L, M, E). Substituting this 
wave function into Eq. (1), one can obtain the radial part equation as follows; 

'iff -^-UvW>-*-m. •*>-(£?!. 

Here fc is the quantum number of the spinor operator K defined by [9, 11] 

K=-P q (Z q -L + l), K^ m = k^ m . (5) 

where S g (= <r q 12x2) and L are the 4-component spin and the orbital angular momentum of the light antiquark, 
respectively. Note that with this quantum number k, one can simultaneously determine both a partial angular 
momentum j q of light degrees of freedom and a total parity P as [11], 

j q = \k\~\, P = jt l (-) lkl+1 - (6) 

It is remarkable that in our approach K can be defined even for a two-body bound system composed of a heavy quark 
and a light antiquark. 

Then, in our model a chiral limit is realized by setting m q = S(r) = 0, in which the corresponding Hamiltonian 
becomes 

HZ hiral = a q -p + V(r). (7) 
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With this chiral Hamiltonian, the radial part equation becomes 

V -dr + l 

dr + l: V 



^chiral ( r ) = E k (chiral) ^chiral ^ (g) 



which means, 



This equation can be converted into the one with — k by a unitary transformation with U = 02, where 02 denotes the 
2nd component of Pauli matrices, 

d V _u - dr - k Au^t ira \r)=E^ chi ^U^t ira \r), (9) 

r r / 

U^t lra \r) = * c _t al (r) and E k Q {chiral) = E ~ k {chiral) . (10) 

That is, the energy is degenerate in ±fc in the case of chiral Hamiltonian, Eq. (7), which corresponds to the degeneracy 
for k = ±1 in the chiral limit. The degenerate hadron mass is given by tuq since Eq.(8) is nothing but an equation 
for a hydrogen atom with mass = and its energy level is proportional to mass which is zero in this case, i.e., 

j^k (chiral) n 

h = U. 

Considering the above observation, mass splitting in our model occurs as follows. 

(1) Start from the chiral limit Hamiltonian, Eq. (7), together with no 1/rriQ corrections with tuq heavy quark mass. 
In this stage, all the masses of 0~, + , 1~ and 1+ states with j q = 1/2 are degenerate; tuq = m(0~) = m(l~) = 
m(0+)=m(l+). 

(2) When the light quark mass m q and a scalar potential S(r), i.e. explicit chiral breaking terms are inserted as 
shown in Eq. (1), degeneracy is partially broken; 77i(0~) = m(l~) and m(0 + ) = m(l + ), which are called heavy 
quark multiplets in [2], because there still remains degeneracy due to the quantum number k [11]. 

(3) Finally by including 1/mQ terms, all the degeneracy is resolved and the mass values in Tables I and II are given. 

Notice that the quantum number k plays an important role in our model to see how the states are classified and 
how the degeneracy is resolved [11]. 

The above procedure may be depicted in Figure 1. 

FIGURE 1: Procedure how the degeneracy is resolved in our potential model. 

k - +1 1+ 



U = 1/2 



0+ 



0, s -> 



/.• I 1 



no 1/mo corrections / 

(m q ^0,S^0) 

{l/rriQ corrections) 

On the other hand, the procedure taken by Bardeen and Hill [3] can be stated as follows. 

(1) Start from the chiral limit of light quark masses but with infinite heavy quark mass instead of finite. In this 
stage, all the masses of 0~, + , 1~, and 1 + states with j q = 1/2 are degenerate just like our model. 

(2) Inclusion of a finite heavy quark mass leads to a partial resolution of degeneracy; to(0~) = m(0 + ) and m(l~) = 
m(l + ), which are called parity doublets. 

(3) Finally by including a finite light quark mass effects and using the modified Goldberger-Treiman relation, all 
the degeneracy is resolved. 

The above procedure may be depicted in Figure 2. 

Comparing these two procedures, one can easily see that our model naturally explains how to resolve degeneracy 
between 0" and 1~ and/or + and 1+ due to the quantum number k and clarify the origin of mass splitting. Namely 
the interpretation by Bardeen, Eichten, and Hill is not the unique way to explain the mass splitting. There still 
remains a way to explain it using a potential model. Moreover, it should be noted that the potential model can give 
not only the mass differences of heavy mesons but also their absolute values, as seen in Tables I and II. 
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FIGURE 2: Procedure how the degeneracy is resolved in [2]. 

1+ 

0+ 

1" 

0" 

(m q + 0) 

IV. CHECK OF THE GOLDBERGER-TREIMAN RELATION 

Before closing this paper, we would like to give some comments on the modified Goldberger-Treiman relation 
claimed by Bardeen, Eichten, and Hill [2]. The modified Goldberger-Treiman relation predicts AM(0) = AM(1), 
where 

AM (0) = M(0+) - M(CT), AM (1) = M(l+) - M(l"). 
Then, one can test it by using the following mass differences for a couple of states. 

(1) D s (cs) states [1] 

AM(0) = £> s ./(0 + , 2317) - D s (0~, 1969) = 348 MeV 
AM(1) = D sJ (l+, 2460) - D* (1", 2112) = 348 MeV 

(2) D(cu(d)) states [7] 

AM(0) = £>(0+, 2308) - D(0~ , 1870) = 438 MeV 
AM(1) = D(l + , 2420) - L>(1~, 2010) = 410 MeV 

(3) K(su(d)) states [10] 

AM(0) = K(0+, 1430) - K(0~, 495) = 935 MeV 
AM (1) = K(l+, 1335) - K(l~, 895) = 440 MeV 

where we have given mass of K(l+) =1335 MeV as a simple average of ifi(1270) and ^(1400). L>(0+, 2308) is from 
the Belle data [7] and others are from the particle data group [10]. As one can see here, the modified Goldberger- 
Treiman relation looks to work amazingly well for D s (cs) states. However, it might be problematic for other cases as 
shown in the following, 

(1) the mass of D(0+) is expected to be 2280 MeV, if AM (0) = AM(1) = 410 MeV holds, which is quite different 
from 2308 MeV observed by Belle [7]. Compared to this expectation, it is remarkable that in our case it is 
estimated to be 2304 MeV, as shown in Table I, being very close to the Belle data. 

(2) furthermore in the case of K states, it seems this relation cannot be applied, though its application to this case 
might not be justified because s quark is not heavy. 

Based on this observation, the modified Goldberger-Treiman relation might be correct only in the case of D{cs) and 
there might be no more meaning for other states. To confirm this, it is very important to observe B(0 + ) and B(l + ) 
as well as B s (0 + ) and B s (l + ), though a collective mass of orbitally excited (L = 1) B meson states was reported a 
couple of years ago [12]. We strongly wish that these particles will be observed soon. 



k = 1/2 



m c — ► oo 



J p =0± 



(m c finite) 
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